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1 Introduction. 



The scalar field theory with a cj)'^ interaction has been widely studied for a variety of problems. 
For instance, it underlies the physics of the anti- ferromagnetic phase transition in statistical 
physics and also is the starting point for the Higgs mechanism of the standard model in particle 
physics. From another viewpoint it has been used as a toy model in four dimensions to examine 
fundamental ideas in quantum field theory. One such example is understanding perturbation 
theory at high orders. (For a review of these points see, for example, Q.) In particular the 
fundamental functions of the renormalization group like the /3-function are known to five loops 
in the MS scheme, |^. In other field theories of interest in particle physics such as gauge 
theories four loop results have only become available in the last few years. For instance the four 
loop /3-functions of QED and QCD, which describes the strong interactions, were computed in 
IQ, Q respectively. 

Such calculations are necessary in QCD primarily because it is asymptotically free and 
therefore it describes the physics of high energy collisions involving quarks. As these deep 
inelastic experiments are becoming increasingly more accurate the theoretical input must be 
refined accordingly. Therefore current perturbative calculations in QCD have been focussed on 
the contribution of the anomalous dimensions and coefficient functions of the twist-2 flavour 
non-singlet and singlet operators that arise in the operator product expansion. These were 
originally computed at one and two loops in ^ as a function of the moment of the operator 
n. Recently the first few moments of the three loop dimensions were determined in fs], ^. 
However, the programme is not yet complete as one requires the explicit n-dependence at three 
loops. This is required in order to compute the DGLAP splitting functions, [l^, which are 
given by the inverse Mellin transform with respect to n. In |^] an approximate fitting procedure 
was developed to obtain these functions. This was achieved by using general information on the 
(expected) form of the function in Bjorken x, which is the variable conjugate to n with < a; 
< 1. Also the information from the explicit exact result for the first few moments was used to 
constrain the fitting parameters. In the absence of a full n-dependent result for the operator 
dimensions this strategy has to suffice for present. Clearly there are several potential problems 
in such an exercise. For example, the first few moments would only be expected to give a fairly 
accurate approximation to the splitting functions in the bulk of the x-range. However, the low 
X region which is presently of physical interest would not be well covered since there the higher 
moments give important contributions. So it would be useful to understand to what extent such 
a fitting procedure can be relied upon. 

Whilst the complete calculation in QCD will be a formidable computational task one could 
at least address this and other issues in a toy model. This is the primary motivation for this 
paper. We will compute the anomalous dimensions of the analogous operators in 0^ theory to 
four loops in the coupling constant expansion in the MS scheme. Although this would appear 
to be one order beyond that currently of interest in QCD, it turns out that due to the form 
of the model there can be no one loop contribution. We view this theory as a toy for several 
reasons. First, the nature of the interaction is simpler than in a non-abelian gauge theory and 
hence there will be substantially fewer diagrams to consider and also no tensorial complication 
in the numerator of the Feynman diagrams which occurs when one has fermion and derivative 
couplings. Consequently although we do not expect the results to be as involved as those which 
will undoubtably occur when the full QCD result is realised, those functions of n, such as finite 
sums in n, which will occur will at least mimic analogous structures. Indeed as will be seen, it 
is these forms which will drive the strategy of the fitting procedure. As we will obtain the full 
4-loop results as a function of n we will be in a position to investigate and gain insight into the 
approximation procedure for the splitting functions. 
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This leads us to our second reason why we will regard our work as a laboratory for testing 
ideas. Clearly unlike QCD cp^ theory in four spacetime dimensions is not asymptotically free 
and therefore it does not sense to think of the high energy limit describing particle collisions. 
However we stress that the aim is a mathematical one to endeavour to understand the relation 
of functions of the operator moment n to their Mellin transform which in (/)^ we will regard as 
being on the same footing as the DGLAP functions. In this situation asymptotic freedom plays 
a minimal role. Finally, our second motivation for this study is in relation to understanding 
the content of the operator product expansion of 0^ in d-dimensions. As noted earlier the 
model relates to the perturbation theory of the Heisenberg model which equivalently, in the 
sense of critical phenomena, can be described by the 0{N) non- linear a model in d = (2 + e) 
dimensions. (See, for example, ||l|.) The operator product expansion in that model has received 



wide attention in recent years both perturbatively and in the expansion, |11|. Therefore the 
computation of the anomalous dimension of a set of operators involving derivatives or gradients 
at high order in the related (j)'^ theory is important for providing at least a cross check on future 
perturbative calculations in this area. 

The paper is organised as follows. In section 2, we review basic features of the renormalization 
of composite operators in quantum field theories and then perform the two and three loop 
computations for the anomalous dimensions of the operators in question. This calculation is 
extended to the fourth order in section 3 where we also discuss the technical details of the 
determination of the divergent part of some of the underlying Feynman integrals. These results 
are then used to deduce the corresponding splitting functions in section 4 where the approximate 
fitting procedure is analysed. Our conclusions are given in section 5 and intermediate results in 
the four loop calculations are given in an appendix. 



2 Preliminaries. 

We begin our discussion by reviewing the background to the computation of anomalous dimen- 
sions of composite operators. We will concentrate on one particular set of operators in this 
paper. 

On = (t){x)d^^ . . . d^^(l){x) - traces (2.1) 

which are twist-2 from the phenomenology point of view and symmetric and traceless. The field 
(/) is the basic field of the (/>^ theory which has the lagrangian 

L = lidc^f - A(^2)2 (2.2) 

where g is the coupling constant which is dimensionless in 4-dimensions. As we will be applying 
standard techniques of renormalization such as dimensional regularization in d = 4 — 2e dimen- 
sions and using the modified minimal subtraction scheme, (MS), we must derive the momentum 
space version of the operator. Its divergent structure is deduced from the divergence structure 
of 2-point Green's functions where On has been inserted. Therefore we write 

/ = / 0^Ji-P)Tn{pU...,J{p) (2.3) 

where Tn{p)^^...^^ is a symmetric traceless tensor. It is more convenient to swamp the Lorentz 
indices by contracting with a constant vector and define the object 

r„(p,A) = r„(p)^,...^,A^^..A'^" (2.4) 
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Given the fact that Tn{p) is traceless one can derive the form of T„(p, A). However, in higher 
order calculations we will exploit the first two terms in the construction of recurrence relations 
for subintegrals. Therefore we record, 

T„(p,A) = (Apr - J^l~/\y A\Apr-^ + 0{{A^)\ApT-^) (2.5) 

In QCD calculations [^, 0, |9| the constant vector is usually taken to be a null vector in which 
case ( |2.5| ) terminates at the first term. 



The anomalous dimension we are aiming to compute is Jnid)- In terms of renormalization 
constants this is defined to be 

Inig) = M^^HZ^^Za) (2.6) 
where M is the mass scale introduced to ensure g remains dimensionless in d-dimensions. The 

1/2 

renormalization constant Z^ corresponds to the renormalization of the field (/) and is known 
at five loops in the more general theory with an internal 0{N) symmetry Q. The lower order 
calculations were performed in In our notation defining 

l{9) = 27<^(5) = M^^lnZi (2.7) 

then, for = 1, 

7(5) = la' - \a' + f/ (2.8) 

The remaining renormalization constant Zq defines the renormalization group function joia) 
via 

7o(<?) = M^-^lnZo (2.9) 

which corresponds to the renormalization of the operator insertion in the one particle irreducible 
2-point Green's function. From the nature of the quartic interaction 70(5) will have the 0{g^) 
term as the first non-zero contribution. In the analogous calculation in QCD we note that the 
twist-2 operators are physical operators and therefore their anomalous dimensions are gauge 
independent. However the renormalization group functions corresponding to 7(5) and 70(g) are 
each gauge dependent. Therefore although there is no gauge symmetry in (j)^ it is 7^(5) which 
we will regard as being on the same footing as the gauge independent dimensions in QCD. 



We now detail the computation of the two and three loop contributions to (2.9). These 
results will play an important role in the four loop calculation where they will enter multiplied 
by vertex counterterms for example. Although strictly it is their e-expansion to 0(e) and 0(1) 
respectively for two and three loops which we will need it is in fact possible to compute the 
relevant integrals exactly as a function of d. The two loop contribution to ( |2.9D arises from the 
first graph of fig. 1 where the line with a dot on it represents the insertion of the operator 
Tn{p,A). There is a flow of momenta p through the external legs of the Green's function. As 
we are interested only in the ultra-violet divergence of the graphs we compute with massless 
propagators which in turn allows us to exploit massless integration techniques such as uniqueness 
1 12, 13]. However for the lower order cases one needs only to apply, for instance, the elementary 



chain rule of |12|. Therefore the exact result for the value the 2-loop graph of fig. 1 is 

1/(1, 1, 2/x - 2)i/o„n(2 - /X, 2, 3/x - 4 + n) (2.10) 

where d = 2/i = 4 — 2e, 

z^(q,/3,7) = 7r^a(Q)a(/3)a(7) , i^mnpia, = TT''am{a)an{P)ap{j) (2.11) 
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a{a) = r{fi — a)/T{a) and = — a + n)/T{a). The p and A dependence is easy to 

reconstruct from the integral dimension. Performing the e-expansion of ( 2.10| ) and including the 
symmetry factor of 1/2 gives the contribution to (|2.9D as 



^ (2.12) 



2n(n + l)e 



The three loop calculation requires the calculation of the second graph of fig. 1. As there 
are two graphs with the same topology we have included both in the same diagram with their 
respective names, El and E2. Again it is straightforward to compute each integral exactly using 



the chain rule of [12| and we find that their values are 

El = (zy(l,l,2/i-2))2i/„o„(2,4-2/i,4/i-6 + n) (2.13) 
E2 = z.(l,l,2/i-2)zy„on(2,l,2//-3 + n)z.„On(5-2;U,l,4/i-6 + n) (2.14) 

However each graph has at least one subgraph divergence arising when either one or both 
loops can be enclosed by a box. The divergence arising from these subgraphs are cancelled by 
multiplying the two loop graphs of fig. 1 by the relevant (vertex) counterterm. In this and 
the four loop case we include in our expressions for the final value of the divergent part of the 
graphs the subtraction of these subgraph divergences. (Further background to this procedure 
is well documented in, for example, ||l^.) However it is appropriate to record that for E2 the 
subtraction is 

a2(l)an(2)a„(3/i-4 + n)/e (2.15) 
Hence the full MS contributions to (^^ at 0{g^) are 

1 (2n2-l) 

El = — + —— ^ (2.16) 

3n(n + l)e^ 3n^(n + l)^e 

E2 = 1 I ( I (4n^ + 2n-l) \ 1 ^^.17) 

6n(n + l)e^ I 3n(?i + 1) 6n^(n + 1)-^ / e 

where the finite sums Si{n) are defined as 5/(n) = J2]=i With these values we deduce that 

. -J—y + (2M± _ + (2.18) 

'^^^ n(n + l)^ \n{n + l) 2n^{n + iy r ^ ' 



3 Four loop calculation. 

We now turn to the renormalization of the operator at the four loop level. There are four 
basic topologies to be considered which are illustrated in fig. 2. As in fig. 1 we have indicated 
the location of the operator insertion by a dot and denoted the corresponding graph by the 
symbol beside each insertion. Of these four topologies sets C and D represent integrals whose 
divergence structure and subgraph subtractions can be computed directly by elementary chain 
integrals involving ( |2.1l| ) and we will not discuss them further. Likewise integral A2 is trivial 
to compute. The remaining graphs required some ingenuity and we detail the techniques for 
several as an aid to the interested reader and a potential method for future similar problems. 

First, we consider the integral ^41 and detail its calculation. It can quickly be reduced to the 
2-loop integral illustrated in the first graph of fig. 3 where the factor arising from the integrations 
is 

a„(l + 3e)(a(l))3 

4e2(l - 2e)(n + 1 - 5e)a(l + e)a„(l + 4e) ^ ' ' 
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In fig. 3 we denote the power of the propagator beside each line if it is not unity and the location 
of the numerator (Ay)" by the symbol (n) and an arrow on the line beside the propagator 
involving y. We have used the coordinate space representation of the integrals so that the 
variables of integration are the location of the vertices and not the flow of momentum around a 



loop, |12|. On dimensional grounds the value of this integral is 




and the aim is to deduce Ai{e,n) to the term linear in e. This is because ( |3.1D is 0(l/e^). The 
first step is to relate Ai(e, n) to an integral with all exponents unity or 0(e) which is achieved 



by the integration by parts rule given in 15]. The result is the set of integrals on the right 
t equation of fig. 3. T 
^e the contribution to A 

a(l)a(l + e)a„(l + 2e) 



side of the first equation of fig. 3. The three terms with a line missing are easy to compute 
exactly and give the contribution to Ai{e,n) of 



6e(n + 1 - 4e)a„(l + 3e) 



a„(l) a(l) 



(n - e)a„(l + 2e) (1 - 3e)a(l + 2e) 



(3.3) 



The remaining two diagrams have similar structure. If we denote by bi{e,n) the value of the 
first graph on the right side of fig. 3 then we can obtain a recurrence relation for bi{e,n) with 
respect to n which can be solved. One way of achieving this is to consider the integral without 
the A contraction and decompose the numerator structure . . . y^^ into a sum of symmetric 
and traceless objects defined in |13]. As such objects are independent the resulting 



integral must be a sum of the set y[^^„,^,.) with r < n where the balance of indices is made up 
by including the appropriate number of r]^u tensors. Their coefficients are functions of n with 
the leading term corresponding to 6i(e,n). This is deduced by eliminating the next to leading 
order term between the equations obtained by first contracting with A^^ . . . A^^ and applying 
the differential operator x^d/d/^^ and second by contracting with, say, 7^^*1-1/^". This results in 
the relation 

6i(e,n) = r„[6i(e,n-l) + ci(e,n-l)] (3.4) 
where ci(e,?i) is defined in the second equation of fig. 3 and 



(n + /i - 2) 
(n + 2^ - 3) 



(3.5) 



Clearly the solution to (^) is 

61(6, n) = 61(6,0) (nr,) + f^ci(6,A;-l) (Hr,) (3.6) 



\i=i / k=i 



The integral 6i(e,0) has been computed explicitly in [|T6| and has the value 

^' ' 3e2(l-3e)2(l-4e)a(l + 3e) ^ ^ 

in our notation. For ci(e, n) the second term is 

1 a(l)a(l + e)an(l) 

3(n-e)(l-3e)(n + l-3e)(n + l-4e) a„(l + 3e) 

whilst the first and third are similar in form where, for example, it is simple to derive the integral 
representation for the first graph of the second equation of fig. 3 as 

3(l-2eKl-3e)ra+t) [ ^"^^^ " ^^"^^ ^ ^^^"^^^^ " ^^^^^^ " ''^^ ^'''^ 



(3.8) 
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Since we are only interested in the divergent part of Al it is straightforward to expand the 
integral in powers of e to the 0(e) term and obtain for it alone 



ds s 
1 



^(1 



\l-2e 



dt t 



l-3e. 



l-t)-'[l-str 



Si{n + 2) 



+ 



(n + 1) 
e 

(n + 1) 



(n + l)(n + 2) 

3S'i(n + 2) 



3 + 5i(n + 2) + 



25^(n + 2) 5S2(n + 2) 



Hence we have 



bi{e,n) 



(n + 2) (n + 2) (n + 2) 

(l-2e)(l-5e) an(l)a^(l + e)a(l) 

3e2(l - 3e)2(l -4e){n + l- 2e) a(l + 3e)a„(l + e) 



+ 0(e2) (3.10) 



+ 



3(n + l) 



2-S2(n + l) 



Si{n + 1) 



(n + 1) 



(3.11) 



Assemblying all the relevant contributions and factors like 
we deduce that the divergent part of Al itself is 



and expanding in powers of e 



1 



+ 



1_ ( Si(n) (Sn^ + lln + l) 

8n(n + l)e3 ' ^ l^4n(n + 1) 8n2(n + 1)2 

1 / (252(n) + 5i2(n)) (5n2 + lln + l)Si{n) 



+ 



4n(n + 1) 



+ 



4n2(n + 1) 



+ 



(19n^ + 68n3 + 85n2 + 12n + 1) 
8n3(n + 1)3 



0(1) 



(3.12) 



Performing the simple operation of subtraction of subgraph divergences yields the result for Al 
which we have recorded in eqn ( [A.ll ). 

Another technique was also used to compute several of the other graphs of fig. 2 and we 
detail this for graph Bl. After performing the two elementary loop integrations we obtain 
a simple two loop graph. Rather than considering this graph we examine the more general 
graph illustrated in fig. 4 where we have introduced an extra parameter a and the (analytic) 
regularization 5. Denoting the value of this graph by Bi{e, o, 5) we record that for the problem 
in hand we require i?i(e, 1,0). The presence of a general a and 5 will allow us more freedom 
in using various calculational techniques such as uniqueness to obtain expressions to determine 
the coefficients in the power series. 

The first step is to apply the operator Af^d/dp^ to the diagram to yield the two graphs Ii 
and I2 of fig. 4 where 

Bi{e,a,6) = [aeli + (1 - ,5)/2] (3.13) 

(3 + a)e 

Whilst the former is finite with respect to e the second is divergent. It is determined by inte- 
grating by parts on the integral with the same exponent structure on each line except that 1 
is subtracted from the top right line and added to the bottom right. Substituting the result of 
this operation into ( 3.13 ) we obtain 

Bi{€,a,S) 



1 



(3 + a)e 



[il-5)J2{e,a,6) + F(e,a,(5)] 



(3.14) 



where J2 is integral given in fig. 4 and F{€, a, 5) represents the sum of the remaining graphs 
which are finite with respect to e. As J2(e,a,5) is a set of elementary chain integrals all that 
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remains is the calculation of F{e, 1,0) whose 0(e) term is needed due to the pole in e in ( 3.14| ). 
This is achieved by expanding in powers of a and e 



Fie,a,6) = fo + (/i + a/2)e + 0{e^) (3.15) 

and then determining the constant coefficients /j. Rewriting F{e,a,S) from ( p. 141 ) as 

F{e,a,6) = {3 + a)eBi{e,a,6) - {I - 6)J2{e,a,S) = D{e,a,6) (3.16) 

then it can be evaluated for special cases of a to give 

F(e,l,0) = 2D{e,0,0) - D{e,-1,0) (3.17) 

In the explicit evaluation of the integrals, however, there is a potential problem when using 
uniqueness to integrate the right triangle of i?i(e, — 1,0). To avoid the appearance of factors 
such as r(0) in both numerator and denominator one needs the temporary regulator 6 and so 
(|3.17| ) ought to be replaced by 

F(e,l,0) = [2D{e,0,5) - D{e,-l,d)]s^o (3.18) 

As an aid we record the results of these intermediate steps as 

a(l)a(l + e)a„(l + e)ai(2) 



J2(e,l,0) 
l?(e,0,0) 



a(2e)a„+i(l + 2e) 
a3(l) / a„(l + 2e) a„(l + e) 



a(e) Van+i(l + 3e) a„+i(l + 2e) 



D(e,-l,0) = + ^ (3,19) 

^ ' ' ^ (l-e)(n + l-3e)2an(l + 2e) ^ ' 



With these expressions it is a simple exercise to verify that the result ( [A.4| ) is finally obtained 
for graph Bl. 

We have used either of these algebraic techniques outlined here to evaluate the remaining 
divergent parts of the integrals. For completeness the results for each are recorded in the 
appendix where the subgraph divergences have been removed. In several of the results the 
additional finite sum Ki{n) = Y^'^=i{~'^)^^^^^ / f occurs. In ( |A.6| ), for instance, it arises in the 
computation of an intermediate two loop integral where the exponents on the lines diagonally 
opposite each other are 2 and 1 — e respectively and the central line has zero exponent but an 
(n) insertion. However, its e-expansion is elementary to compute using the rule for chains after 
the application of the binomial expansion. Also in the graph A2> after completing the first two 
elementary loop integrals one is left with a two loop integral multiplied by 1/e. This integral, 
which is finite with respect to e, has been evaluated explicitly in |]l^] and involves K2{n). 

Finally it is a simple exercise to assemble all the contributions to the four loop result and 
allowing for the symmetry factors of the graphs which are given in table 1, we find the 0{g'^) 
term of ( |2.18| ) is 

352(n) ^ Sfin) (89n2 + 53n - 18)^1 (n) 



n(n + 1) n2(n + 1)2 6n2(n+l)2 
_ (n^ + n - 4)7^2 (n) (265n^ + 28071^ - SGn^ - 39n + 33) 

n2(n+l)2 + 12n3(n+l)3 ^ ^ 
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We summarize the results of this and the previous section by recording our 4-loop MS value for 

Inig) as 

(n-2)(n + 3) 2 f^SUn) (n^ + 2n^-39n'^-16n + 12) 



6n(n + 1) \n{n + l) S'n?{n+1)'^ 

3S'2(n) ^ 2Sf{n) (SDn^ + 53n - 18)5i(n) 



\n{n + l) n{n + 1) 6n2(n + 1)^ 

(n^ + n -4)K2(n) {265n^ + 280n^ - 36n^ - 39n + 33) 65\4 

n2(n+l)2 + 12n3(n+l)3 Slj^ ^ ^ ^ 

One non-trivial check on this result is that it vanishes for n = 2. In that case the original operator 
corresponds to the energy momentum tensor which is conserved in the quantum theory. It is 
well established that the anomalous dimension of non-anomalous conserved physical currents 
vanish to all orders in perturbation theory. (See, for example, [T^] .) 



4 Mellin transform. 



We now apply the results of the previous sections to the problem outlined in the introduction. 
Again we emphasise that the aim of the exercise is to investigate the mathematics underlying 
the fitting to the Mellin transform of the anomalous dimensions given only knowledge of the first 
few moments. In QCD it is this transform which corresponds to the DGLAP splitting functions, 
|llO|| . These are a measure of the probability that a parton decomposes into other partons ie 
quarks and gluons. The variable conjugate to the moment n is Bjorken x which represents the 
momentum fraction carried by the parton in the nucleon. Therefore it is restricted to lie in the 
unit interval. Though in general the domain of x is the half-line. To approximate the three loop 
splitting functions in QCD the authors of |^] took a set of trial functions of x whose properties 
were consistent with expected general behaviour of the splitting functions in the unit interval 
and fixed the unknown coefficients using the available exact 3-loop moments. For example, for 
the non-singlet case these expressions are known for n = 2, 4, 6, 8 and 10. Of course such 
an exercise is unnecessary if the full n-dependent result was available. In |^] as a prelude to 
exploring the third order correction the two loop result was examined. That had the advantage 
of knowledge of the exact result of the operator dimension for all n and hence the splitting 
functions for all x. One disadvantage is that the second order result has a simpler structure of 
functions of x than those which would arise at next order. Also in this second order case the 
result is unlikely to be unambiguous as only a few moments are sufficient to isolate which of the 
set of trial functions play the important role in, say, the small x behaviour. 

Having summarized the status of the QCD calculation we can use our results at four loops 
in (j)'^ to try and gain a more detailed insight into improvements to the fitting. First, some 
elementary remarks are in order. It is important to note that we are not claiming there is a 
parton interpretation in (f)'^ theory. The lack of asymptotic freedom counts against this. Second, 
although the QCD calculations are at third order and have computed 7n(5) to four loops, it is 
in effect a third order calculation since there is no one loop contribution. Therefore it is this 
fourth order result which we will regard as being on the same footing as the three loop QCD 
result although we are not claiming that it has as complex a form. We begin by defining the 
Mellin transform of a function /(x) as 

M[fix)] = dxx'^-^f{x) (4.1) 
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where we restrict the function to exist on the unit interval. In QCD the operator dimension 
is gauge independent though the field dimension and the renormalization of the bare operator 
are gauge dependent. As there is no gauge symmetry in (j)'^ theory rather than take the full 
dimension 7n(fl') we concentrate on the part 70(5). This is justified by the fact that 7(5) is 
n-independent and therefore its contribution to the x-dependence will be purely 'y{g)5{l — x) 
which is uninteresting for our study. With the definition (^]^) it is straightforward to verify that 
the splitting functions are given by 

00 

P{g,x) = Y.Pn{x)g^ (4.2) 



where 



P2{x) = (4.3) 

^3(2;) = -2xln(x) - 2(1 - x)ln(l - x) - i(l - x)(3 ln(x) + 10) (4.4) 

Pi{x) = 8(1 -x)Li3(-x) - 4(1 -x)ln(x)Li2(-x) - 3(1 - x)V'"(l) 

2 

- 4xln(x)V''(l) H xln^(x) + 9(1 + x) ln(x) ln(l + x) 

3 



+ 9(1 + x)Li2(-x) - 2(l-x)Li2(x) - 4xLi2(x) 

6 



89 

+ (1 -x)ln(x)ln(l -x) + 2(1 -x)ln2(l -x) + — (l-x)ln(l 



93 

+ 13xV''(l) + 3(l-x)^'(l) + — (1-x) 

23 11 27 9 

H ln(x) + — In^(x) — — xln^(x) + — xln(x) (4.5) 

2 8 8 2 

In these expressions Li„(x) is the polylogarithm function whose properties are well known, [p!9|], 
and ^{x) is the logarithmic derivative of the Euler T-function. To gain some idea into the form 
of these functions we have plotted them in figs 5-7. 

As a preliminary to the four loop analysis we first summarize the method for the three loop 
case since the two loop example is clearly trivial. The basic idea is to postulate a basis set of 
trial functions of x with arbitrary coefficients. These are fixed by using (|4.1| ) to obtain a function 
of n which can be evaluated for the first few moments and compared with the analogous value 



of ( 2.18 ). Once the coefficients have been fixed one can compare the approximate function with 
the exact result to see how well it covers the function in the unit interval. It is worth pointing 
out that the trial set of functions ought to contain a few which have logarithmic dependence 
given the nature of the exact result and its form in fig. 6. Therefore some guidance on how 
good the fit is, is that it should be better than taking a trial set which gives a polyniomial in x 
of the appropriate degree. Given this criterion we discovered that the set 

{ 1, X, ln(x), ln(l - x), X ln(l - x) } (4.6) 

gave a very good approximation. We have plotted this approximation with the exact result 
in fig. 8 where the straight line is the exact result and the dots represent the approximation. 
Interestingly the discrepancy at x = is about 0.12. Clearly to achieve this accuracy we have 
used only five moments. However, comparing the functions of ( [4. 6)) with the exact function one 
observes that xlnx is absent. This illustrates a potential pitfall in fitting approximations to 
exact results. One feature of the three loop case to emerge was the necessity of having functions 
in the trial set which reproduced some of the functions of x which appeared in the exact result. 
We repeated this exercise with a variety of sets either with fewer functions and therefore fewer 
moments or replaced several of the entries in (|4.6| ) with other functions. However the result of 
fig. 8 represents the best fit. 
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We have repeated this analysis for the 4- loop case. From the form of the exact result it turns 
out to be a harder exercise especially in covering the small x region well. Again with the criterion 
that a fit must be better than a polynomial, we managed to gain reasonable approximations 
with two trial sets. To achieve this we needed to extend the set to seven elements and include 
higher powers of the logarithm. Even with this, however, an answer as accurate as three loops 
could not be achieved. These sets were 

{ 1, X, ln(x), In^(x), Li2(x), ln(2;)Li2(2;), X ln(x) } (4.7) 

{ 1, X, ln(x), ln2(x), Li2(x), ln(j;)Li2(x), Li2(-x) } (4.8) 

and we have plotted each beside the exact results respectively in figs 9 and 10. Although set ( |4.8| ) 
appears to deviate in the middle x range it is clearly more accurate for small x than ([4.7| ). This 
is due to the alternating form of the expansion of the dilogarithm for small x. The respective 
discrepancies at x = from the exact and approximate answers were 2.5 and 0.5. 



5 Conclusions. 



We have constructed the MS 4-loop anomalous dimensions for a set of twist-2 operators which 
are analogous to operators which occur in the operator product expansion used in deep inelastic 
scattering in QCD. In the course of examining the inverse Mellin transform for the four loop 
case we found that it was important to isolate those functions of x and include them in the fit 
that reproduced similar forms to the basic finite sums that appear in the exact answer. In the 
absence of a full exact answer for QCD it would therefore seem to us that one could endeavour 
to improve the fits of the anomalous dimensions, in the low x-region by trying to deduce the 
type of finite sums which would occur at three loops. An insight into this could be obtained by 
generalizing some of the series that already occur in the exact two loop results. For instance, 
the sums Si{n) occur at Z-loops for / = 1 and 2 and it would appear odd if 53(n) did not occur 
in the 3-loop result. Therefore another structure to study might arise from generalizing those 
finite sums which only occur in the two loop result. The other feature which we observed in our 
analysis was that to maintain the accuracy one would require several more exact moments for 
the anomalous dimensions. 
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A Values for 4- loop integrals. 



In this section we record the values of the divergent parts of the integrals which arise in the 
computation of the 4-loop anomalous dimension in MS. The label refers to the diagrams of fig. 
2 though the values given here correspond to those with the subgraph divergences removed. We 
found 



Al = 1 , f Sr{n) _ (7n^ + 5n-l) \ J_ 

24n(n + l)e3 \l2n{n + l) 24n2 {n + lf ) 

^ ( S2{n) Slin) (7n2 + 5n- l)5i(n) 



6n(n + l) 12n(n + l) 12n2(n + l)2 

(25ra^ + 36n^ + 7n^ - 4ra + 1) \ 1 
^ 24n3(n + l)3 ) ~e ^ ' 

^2 _ 1 (5ra=^ + 3n - 1) (9n^ + Sn^ - 3n^ - 2n + 1) 

12n(n + l)e3 ~ 12n2(n + l)2e2 I2n^{n + lfe ^ ' ^ 

^3 = — 57 ^ i\o (■^■3) 



^ 1 . / 5i(n) (7n2 + 5n-l) \ 1 

24n(n + l)e3 \^12n(n + l) 24n2(n + l)2 J 

/ S2{n) Sf{n) _ {7n^ + 5n-l)Si{n) 

^ l^3n(n + l) ^ 12n(n+l) 12n2(n + l)2 

(25n^ + 3G»3 + 7»^-4» + l) \ 1 

24n3(n + l)3 ^ ^ ^ 

iJ2 = 1 I 1 (lln2 + 5n-3) \ 

8n(n + l)e3 e2i^l2n(n + l) 24n'^{n + lf J 

( S2(n) S\(n) _ gi(n)(n2-n-l) 

i^l2n(n + 1) ~ 6n(n+ 1) 12n2(n+ 1)2 

(13ra^ + 4n3 - 11^2 - 2ra + 3) \ 1 
^ 24n3(n + l)3 J e ^ ^ ^ 

1 / 5i(n) (4n2 + 2n- 1)\ 1 



53 = 



+ 



12n(n + l)e3 \^6n(n + l) 12n2(n + l)2 j 

^ ( S2{n) ^ Slin) 5i(n)(4n2 + 2n-l) 
l6n(n + l) 6n(n + 1) 6n2(n + l)2 

K2(n) ^ (4n4-5n2-n + l) \ 1 



2n(n + l) 12n3(n + l)3 



^ (2n^ - 1) _ (4n3 + 3n2-n-l) 

4n(n + l)e3 4n2(n + l)2e2 4n3(n + l)3e ^ ' 

ro = r-^ = 1 , / Sijn) (4rz2 + 2n-l) \ 1 

12n(n + l)e3 ^ Un(n + 1) 12n2(n + l)2 i e2 
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52 (n) ^ S^n) _ 5i(n)(4n^ + 2n-l) 
6n(n + l) 6n(n + 1) 6n2(n + l)2 

(4n'^-5n2-n + l)\ 1 
12n3(n + 1)3 J e 



1 (13n2 + 9n - 2) 

16n(ri + l)e2 32n2(n + l)2e 

1 / 51(71) (3n + 2)(3n + 1)\ 
16n(n + l)e2 V8n(n + 1) 32n2(n + l)2 ) 
1 (n^ - n - 1) 
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Graph 


Al 


A2 


A3 


Bl 


B2 


B3 


CI 


C2 


C3 


Dl 


D2 


D3 


Factor 


1 


1/4 


1/2 


1 


1/2 


1/4 


1/8 


1/2 


1/4 


1/6 


1/6 


1/4 



Table 1. Symmetry factors for each graph defined in fig. 2. 
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Figure Captions. 



Fig. 


1. 


Two and three loop graphs contributing to 'Jnid)- 


Fig. 


2. 
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tour loop graphs contributing to 7n(5j- 


Fig. 


3. 
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Intermediate steps in the calculation ot Al. 


Fig. 


4. 


Intermediate integrals in the calculation ot Bl. 


Fig. 


5. 


Two loop splitting function. 


Fig. 


6. 
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ihree loop splitting function. 


Fig. 


7. 


Four loop splitting function. 


Fig. 


8. 


Approximation to the three loop splitting function. 


Fig. 


9. 


Approximation to the four loop splitting function with first trial set. 


Fig. 


10. Approximation to the four loop splitting function with second trial set 
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ci(e,n) = 




Fig. 3. Intermediate steps in the calculation of Al. 
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Bi{e,a,6) = 



hie, a, 5) = 




Fig. 4. Intermediate integrals in the calculation of Bl. 
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Fig. 6. Three loop splitting function. 




Fig. 8. Approximation to the three loop splitting function. 
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Fig. 9. Approximation to the four loop splitting function with first trial set. 
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Fig. 10. Approximation to the four loop splitting function with second trial set. 



